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April 1, 2013 
Abstract 

The dynamics of microphase separation and the orientation of lamellae in diblock copolymers is 
investigated in terms of a mean-field model. The formation of lamellar stmctures and their stable 
states are explored and it is shown that lamellae are stable not only for the period of the structure 
corresponding to the minimum of the free energy. The range of wavelengths of stable lamellae is 
determined by a functional approach, introduced with this work, which is in agreement with the results 
of a linear stability analysis. The effects of the interaction of block copolymers with confining plane 
boundaries on the lamellae orientation are studied by an extensive analysis of the free energy. By 
changing the surface property at one boundary, a transition from a preferentially perpendicular to a 
parallel lamellar orientation with respect to the boundaries is found, which is rather independent of 
the distance between the boundaiies. Computer simulations reveal, that the time scale of the lamellar 
orientational order dynamics, which is quantitatively characterized in terms of an orientational order 
parameter and the structure factor, depends significantly on the properties of the confining boundaries 
as well as on the quench depth. 



1 Introduction 

During microphase separation diblock copolymers 
can form various nanoscopic structures like lamel- 
lae, cylinders, spheres or bicontinuous gyroids, 
depending on their composition. These self- 
organized periodic nanoscale patterns in bulk ma- 
terials as well as in block copolymer (BCP) films 
attract great attention because of interesting phe- 
nomena in these systems and promising applica- 
tions in nanofabrication, see, e.g., the reviews.-"- 
For lamellar structures the typical width of 
lamellae is of the order of the length of two cova- 
lently bonded polymers and range from 10 ■ ■ ■ 100 
nm. They are readily tunable by varying the 
molecular weights of both blocks of the copoly- 
mer. In bulk materials lamellae are locally ordered 
but on larger length scales one finds a random ori- 



entational order.— Near a substrate lamellae are 
oriented parallel or perpendicularly onto it. The 
orientation is a direct result of a surface and inter- 
facial energy minimization. When coating a sub- 
strate with the same material as one block of a di- 
block copolymer, this block is selected and lamel- 
lae orient parallel to the substrate. Covering the 
substrate with a thin film of a equimolar random 
copolymer, its interaction with the two dissimi- 
lar BCP blocks can be balanced, so that the sub- 
strate behaves neutrally.-^*^ In this case the lamel- 
lae orient perpendicularly near the substrate but 
the lamellae orientation in the plane of the sub- 
strate is disordered at large length scales. 

Various strategies are investigated to achieve 
a long-range orientational order of the lamel- 
lae in BCP films by the application of electric 
fields,—"— shear flow,—"— directional solidifica- 
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tion,-22r^ use of topographically-^^— or chemi- 
cally ^'^^ patterned substrates. 

By chemical patterning of substrates with a pe- 
riodicity close to the lamellae width, a long range 
orientational order of lamellae can be induced.— 
Such long range order can be achieved even in the 
case of small mismatches between both periodic- 
ities. This raises the question, whether straight 
lamellae are also stable at a wavelength apart from 
the optimal one at the minimum of the free energy, 
Ag. Does for lamellae in BCPs also exist a con- 
tinuous wave number band around 2n/Xe, similar 
as in other common pattern forming systems,-^^*^^ 
even further below the critical temperature of mi- 
crophase separation? This question is investigated 
in Sec. [31 where we introduce a general method 
for the determination of stable wave number bands 
in pattern forming systems with a potential, like 
BCPs. We determine by this method the stable 
wave number range of lamellae and we find also 
perfect agreement with the results of a standard 
linear stability analysis of straight lamellae (see 
also Appendix lAl and iBl) . In Sec. [3] we calculate in 
addition analytical solutions for the lamellar struc- 
ture in the weak as well as in the strong segregation 
limit and present also analytical results for the sta- 
bility range in the weak segregation limit in part 

m 

To induce a long range orientational lamellar or- 
der in the plane of BCP films with the lamellae per- 
pendicular to the substrate, the film may be con- 
fined between two lateral boundaries at small and 
medium distances as in Refs.— >^ It depends again 
on the surface preparation of the lateral bound- 
aries, whether lamellae become oriented by en- 
ergetic reasons either parallel or perpendicularly 
with respect to them. The homogeneous lamellar 
structures in such quasi two-dimensional systems 
confined between two boundaries are analysed in 
Sec. S] and Appendix O We determine for vari- 
ous selectivities at the boundaries the lamellar ori- 
entation corresponding to the lowest free energy. 
For nonsymmetric boundary conditions we find in 
two dimensions as a function of the surface selec- 
tivity of one boundary a transition from a prefer- 
entially perpendicular to a parallel orientation to 
the boundaries, which is rather independent of the 
distance between the boundaries. The dynamical 
evolution, including the coarsening and the devel- 



opment of the orientational order of lamellae be- 
tween two boundaries is investigated in Sec. |5] 
The final Sec. [6] includes besides a discussion also 
experimental suggestions. 

2 Model equation 

Microphase separation in an incompressible AB- 
diblock copolymer melt is described in terms of 
a time-dependent Ginzburg-Landau model for the 
conserved mean-field order parameter ^f{r,t) ~ 
<^A(r, — (^B(r,?), with the local concentrations 
(^A,B of the components A and B. 

Spatial variations of the order parameter involve 
a spatial dependence of the chemical potential 
ll{r,t) and the mass current j(r, which deter- 
mine the dynamics of ^f{r^t) via the continuity 
equation 



dt 



-Vj(r,0. 



(1) 



Gradients of the chemical potential drive the mass 
current 

j(r,0 = -MVAi(r,0, (2) 

with the Onsager coefficient M(> 0)-^^ that de- 
scribes the mobility of the monomers A with re- 
spect to B. The functional derivative of a free en- 
ergy functional F{ i//} with respect to the order pa- 
rameter determines the chemical potential 



(3) 



and via Eq. (O and Eq. Q also the dynamics of 



dt 



MV- 



(4) 



The free energy acts as a global Lyapunov 

functional and is always decreasing with time to- 
wards its minimum: 



dt 



dF dy/ 
5xir dt 



dr = M 



-M 



.5F_ 



dr<0. 



— V2— 1 dr 



(5) 



Employing a generalized random phase approx- 
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imation, the bulk free energy functional for di- 
block copolymers was derived by Leibler^? in the 
weak segregation limit, which is applicable to the 
slightly quenched regime of microphase separa- 
tion. Here we use the extended free energy intro- 
duced by Ohta and Kawasaki,— 



K 



dr 



D 



V V 

which also includes the strong segregation limit, 
applicable to the deeply quenched regime. The 
functional comprises the temperature independent 
phenomenological constants u, K, D and \[r = 
(i//) is the spatial average of the order parameter 
V/(r, t). bis the temperature dependent control pa- 
rameter of the model and for b > be microphase 
separation sets in. The second term in Eq. (|6]) with 
the double integral covers the long-range interac- 
tion due to the connectivity of the subchains and 
the Green's function G(r, r') satisfies Laplace's 
equation V^G(r,r') = —d{r — r'). 

The coefficients of the mean field free energy in 
Eq. © can be related to microscopic models un- 
der the following assumptions: All chains have the 
same index of polymerization A'^, are composed of 
by the same number Na{Nb) of monomers of type 
A (B), with = Na+Nb, and have therefore the 
same composition f = Na/N. Both blocks have 
the same Kuhn statistical segment length /. 

The control parameter b in Eq. ^ is related to 
the Flory-Huggins interaction parameter x as fol- 
lows,— 



b = 2x 



2NfHi-fy 



(7) 



where s{f) is of order unity and depends on ap- 
proximations. -i^^ Microphase separation occurs 
for X > %c {T < Tc) and the temperature depen- 
dence of the Flory-Huggins interaction parameter 
X is taken as 



X=A + B/T 



(8) 



where the coefficients A and B are determined 
from experiments.— The parameter m > in 



Eq. Q remains undetermined in the strong seg- 
regation limit.— For the lamellar structure in the 
weak segregation limit u can be identified with the 
vertex function u = r4(0,0), which depends on 
the composition / and the polymerization index 
A/^.iiiM The parameter ^ > describes the inter- 
facial energy between A and B domains and it de- 
pends on the composition / as follows:— 



+- / / G(r, r') [Y{r, t) - y] [v/(r', t) - xjf]drdr', (6) 



K 



4/(1-/) 



(9) 



The parameter 



D 



/W(i-/)2 



(10) 



is positive and decays with the polymerization de- 
gree A'^.— Similar functional dependencies of K 
and D have been obtained in Refs. ^^-^^ except the 
difference in numerical factors due to the use of 
different models for the polymer chain (and there- 
fore different expressions for the radius of gyra- 
tion). Besides the derivation of the interaction pa- 
rameters from microscopic models it is also possi- 
ble to determine them by fitting data from scat- 
tering experiments obtained immediately after a 
quench (see, e.g.,— "^"^ and references therein). 

With the functional given by Eq. ^ and Eq. dH) 
the nonlinear evolution equation of the order pa- 
rameter V^(r,?) follows: 



MkeT [V^{-b\j/ + u\i/^ -KV^Xj/) 

-D{xif-n- (11) 



With the length scale = \fK and the time scale 
T = (^^/ (MhsT) one may introduce with r = i^r' 
and t = xt' the dimensionless variables r' and t' . 
Using in addition the rescaled order parameter 
i//' = -\/mV/, one obtains the dimensionless form of 
the equation (primes are omitted) 

dtXif = W^{-£Xif+\if^-W^\if)-a{\if-p), (12) 

with the dimensionless parameters: 

£ = b, a=DK, p = ^/U^f. (13) 

According to Eq. and Eq. (flOl) one has the scal- 
ing a oc and the limiting case a = corre- 
sponds to the Cahn-Hilliard equation— describing. 
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e.g., phase separation in polymer blends. 

The bulk free energy in dimensionless form is 
given by 



Fb{¥} 



oc. 



+-Mr)(v/-i8) 



dr 



(14) 



where we have introduced for practical reasons the 
auxiliary function 

h{r) = |G(r,r')[v/(rO-/3]Jr'. (15) 



The function h{r) fulfills Poisson's equation 
v2/^(r) = -(v/(r)-/3), 



(16) 



with the boundary condition along the direction n 
normal to the surface of the volume V, 



n-V/z(r) =0 



(17) 



that follows from the conservation of the order pa- 
rameter: /(i//-/3)Jr = 0. 

2.1 Effects of boundaries 

To study unconfined systems periodic boundary 
conditions of the order parameter i// can be ap- 
plied in each spatial direction. For confined sys- 
tems one needs for the fourth order Eq. (1T2| | two 
boundary conditions. The first boundary condition 
follows from the conservation of the order param- 
eter / ( I// — /3 ) Jr = corresponding to a zero flux 
at the boundary: 



n- V (-£!//+ v/^-V^i//) =0 



(18) 



The free energy in a confined system includes a 
surface contribution F^, which can be written in 
dimensionless form, 

P^ = \j s{¥-WsfdS, (19) 



with two phenomenological parameters g and 1//$. 
^ > is a measure of the strength of the interac- 
tion of the block copolymer with the surface and 



V/5 is the preferred difference between the con- 
centrations of A and B at the surface. The case 
g = const, and y/s = const, corresponds to a ho- 
mogeneous surface and non-constant g = g{S) and 
V/(5) model patterned surfaces as e.g. described in 
Ref.— The second boundary condition is derived 
from the local equilibrium condition of the total 
free energy F = F}-, + F^ at the surface 



dF_ 



Ai + [n-Vv/ + ^(V/-V/,)]5- (20) 



Since the bulk is in equilibrium with the surface 
one has the second boundary condition: 



n-Vv/ + g(v/-v/5) =0. 



(21) 



Note that the surface energy in Eq. (1191 ) is equiv- 
alent to the expression proposed in Refs.,— 



j (^-H.xif+^xir^^dS, (22) 



where the "field" Hi is related to the difference of 
chemical potential between A- and B blocks at the 
surface and the parameter ai is related to the so- 
called "extrapolation length", ai ~ that de- 
scribes the ability of the surface to modify the lo- 
cal interaction parameter X-"^^ In our notation one 
has Hi = gyrs and ai = g. The usual situation for 
diblock copolymers is the so-called "ordinary tran- 
sition" with 5 >0 and therefore ai > (or g > 0). 
The surface modifies the local monomer interac- 
tions only within a thin surface layer of thick- 
ness i^(/) with / as the Kuhn statistical segment 
length. In this situation the local interaction pa- 
rameter X at the surface is smaller than in the bulk 
and there is no ordering transition in the range 
T > Tc for Hi = 0. For Hi ^0 finite values of the 
order parameter i/A 7^ are already induced beyond 
the critical temperature, T > Tc. Walls with the 
property Hi^O axe so-called selective boundaries 
and with Hi =0 so-called neutral boundaries. 

Although we restrict our investigations to the 
most relevant ordinary transition we shortly men- 
tion for completeness also another type of transi- 
tion. The so-called "surface transition" occurs for 
ai < 0, which corresponds to 5 < 0. In this case 
the local interaction parameter x is greater than in 
the bulk and even for //i = an ordering transition 
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may be observed at T > Tc near the surface. 

3 Unconfined system: Periodic 
solutions 

Above the onset of microphase separation a per- 
fect lamellar order of block copolymers is de- 
scribed by periodic solutions of Eq. (fT2]) and their 
properties in unconfined systems are investigated 
in this section for symmetric diblock copolymers, 
i.e. /3 = 0. An analytical approximation of the 
amplitude of the spatially periodic solution imme- 
diately above onset is given in Sec. 13. 2[ A method 
for the determination of the stability boundaries of 
periodic solutions, which works close to and even 
far beyond onset of microphase separation is pre- 
sented in Sec. l3.4[ It is based on an analysis of the 
free energy functional evaluated for periodic solu- 
tions. A conventional linear stability analysis of 
nonlinear periodic solutions, cf.,— is presented for 
the present system in Appendix lA] and iBl 

3.1 Onset of microphase separation 

The homogeneous phase of a symmetric diblock 
copolymer melt is described by a vanishing order 
parameter: i// = 0. This basic state becomes unsta- 
ble with respect to small perturbations ^f{r,t) ~ 
gfTfgik r ^j^ej^ tj^g control parameter £ is raised 
beyond its critical value £c, corresponding to a 
quench below the critical temperature Tc of the 
diblock copolymer melt. In this case the growth 
rate o of the perturbations becomes positive and 
microphase separation sets in. Since Eq. (fT2] ) is 
isotropic in space, o depends only on the modulus 
of the wave vector, k = |k|, which is determined 
by the linear part of Eq. (fT2]) : 

o{k)^k^{£-k^)-a . (23) 

The neutral stability condition <j{k) = yields the 
neutral curve 

£f^{k)=k^ + a/k^ (24) 

and the basic state, i// = 0, is unstable above eN{k) 
with a{k) > 0. The critical value of the con- 
trol parameter, £c, and the critical wave number. 



kc, are both obtained from the extremal condi- 
tion deN{k)/dk = at the minimum of the neutral 
curve: 

£, = 2v^, kc = a^^^. (25) 

The wave numbers along the left and right part of 
the neutral curve are given as a function of the con- 
trol parameter £ by the expression 

kl{£) = \(£T^e'-e^^ (26) 

and the wave number k,„ at the maximum of the 
growth rate o{k) increases with £ > £c as follows: 




For further discussions also the reduced control 
parameter r = £/£f — 1 and the reduced wave num- 
ber k = k/kc are useful. Then the rescaled neutral 
curve takes the following form 

r,{k) = (28) 

and the critical values of the control parameter and 
the wave number are given by 

r, = 0, h = \. (29) 

The reduced wave numbers along ryv(^) are 

kl = \+rT^r^ + 2r , (30) 

and the reduced wave number at the maximum of 
the growth rate oijk) is k^ = ^/l+r. 

Microphase separation and spatially periodic so- 
lutions of Eq. (fT2l) develop in the range r > and 
£ > £c, respectively. 

3.2 Amplitude equation 

The basic Eq. (fT2)) takes in terms of the reduced 
control parameter r = £/£c — 1 the following form: 

dtW = V2[-2^2^r + 1) V/+ _ v2^] _ ^4^ (3 J) 

Equation (|3TI) is rotationally invariant and there- 
fore the wave vector of a periodic solution can be 
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chosen parallel to the x-axis with k = {k, 0). 

In the range of small values of r > the neu- 
tral curve is still narrow around kc and nonlinear 
periodic solutions exist only for k rather close to 
kc. Small deviations of the wave vector k from 
the critical one, {kc,0), and therefore long wave- 
length (slow) modulations of the periodic solution 
oc exp(?kf ■ r) can be taken into account by a spa- 
tially dependent amplitude (envelope), 



\j/ = A{x,y,t)e'''''' + C.C. 



(32) 



with A{x,y,t) slowly varying on the scale 2n/kc. 
Such a separation into a slowly varying ampli- 
tude and a fast varying periodic part is successfully 
used in a broad class of pattern forming systems, 
as described for instance in Refs.— i^"^"^ . By 
this separation into short and long length scales 
near threshold, r > 0, a further reduction of the 
basic equation (|3T]) to a universal equation for the 
envelope A (x, J, r) is possible^'^'^'^ and allows 
in the weak segregation regime further analytical 
progress, as described in the following. 

The partial differential equation describing the 
dynamics of A{x,yj), the well-known Newell- 
Whitehead-Segel amplitude equation - for the en- 
velope of periodic solutions in isotropic systems, 
can be derived by a multiple scale analysis,' 



ZodtA 



rA + ^^{d^-—d^ 



A-go\A\^A 



(33) 



"^'^ ^°=2^' ^^"^^ 
Eq. (|33]) has in the range r > periodic solutions 

A=Aoe'^Q'+^y^ (35) 

in terms of the deviations Q = k — k^ and P « k^ 
from the critical wave vector {kc,0). 

The rotational invariance of the system allows 
to choose for stationary solutions P — 0, i.e. the 
stationary amplitude Aq can be expressed in terms 
of 2 or ^ = k/kc'. 



2kf. 



The amplitude Aq vanishes along the neutral curve 



rN = 2Qyk^.=2{k-iy 



(37) 



which is equivalent to Eq. ( |28] ) for r > and 
121 « 1- 

The amplitude equation (|33l) has also non- 
periodic, inhomogeneous solutions, as for instance 
described in Refs.— '^'^^^'^ . A{x,y,t) varies in 
these cases on length scales along the x and the 
y direction, which are larger than In/kc. These 
length scales t,\ along the x direction and t,2 along 
the y direction are: 



1 



1/2 



1 



J_N 1/4 
2r 



(38) 



Accordingly, the envelope A (x,};,?) varies perpen- 
dicular to the lamellae (along the x direction) on a 
different length scale than parallel to the lamellae 
(along the y direction), when for instance the en- 
velope decays from its bulk value A <=< ^/r beyond 
threshold (r > 0) to a small value at the bound- 
ary. The ratio between the two length scales is 
^2/ = (r/S)^/"^ and therefore in the weak segre- 
gation regime r > the length ^2 is always consid- 
erably smaller than the length ^i. This difference 
has a strong influence on the orientation of lamel- 
lae near boundaries as discussed in Sec. |4]. 

3.3 Nonlinear solutions 

In extended systems with periodic boundary con- 
ditions spatially periodic solutions of the nonlinear 
equation (fT2l) of wave number k can be represented 
by a Fourier series. 



M 



j=-M 



(39) 



where the coefficients of this series are determined 
numerically, as described in Appendix |A] For a 
truncated ansatz with one mode. 



Xff = aocos{kx) , 
one obtains for the amplitude gq: 



(40) 



ao = ±2\l^[r-rN{k)] 



(41) 
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which becomes in the range \Q\ << 1 identical to 
the expression in Eq. (|36l ). Again uq only exists be- 
yond the neutral curve r > r^^Qc) [resp. £ > eN{k)]. 

3.4 Wave number bands of stable pe- 
riodic solutions 

Spatially periodic solutions in extended pattern 
forming systems are stable only in a subrange of 
the wave number band beyond a neutral curve as 
given for example by Eq. (|24)) .— 

Stationary, spatially periodic solutions may be 
destabilized, for instance, by small perturbations 
with a wave vector parallel to that of the nonlin- 
ear periodic pattern, if the so-called Eckhaus sta- 
bility boundary is crossed. Or, a periodic solu- 
tion may be destabilized by perturbations with the 
wave vector perpendicular to that of the pattern 
(zig-zag instability), or by a combination of both 
types of destabilizing modes (skewed varicose).^ 
Such stability boundaries are determined by the 
condition, that the growth rate of small perturba- 
tions with respect to nonlinear periodic solutions 
vanishes, as described in more detail in Appendix 
lAlandlll 

In systems where the dynamic equation of the 
field XjT can be derived from a functional F{\j/}, the 
Eckhaus stability boundary and the zig-zag stabil- 
ity boundary can be determined by an analysis of 
the functional F in terms of the periodic solutions 
^f{r) . The idea of this method was indicated ear- 
lier— and it is described below. 

3.4.1 Eckhaus stability boundary 

By crossing the Eckhaus boundary, nonlinear pe- 
riodic solutions become unstable with respect to 
small longitudinal perturbations with a wave vec- 
tor parallel to k of the unperturbed pattern.— !^>^ 
The wave vector k = (fc, 0) is chosen along the x- 
axis and in order to determine the Eckhaus bound- 
ary we calculate the free energy F for a slightly 
perturbed solution, which has a slightly "com- 
pressed" periodicity of wave number k + Ak^ and 
free energy F{k + Akx) /2 in one half of the system 
and a slightly "dilated" periodicity of wave num- 
ber k — Akx and free energy F{k — Akj^)/2 in the 
other half. The free energy of the perturbed peri- 



odic solution is then as follows: 

F{k) = ^-[F{k + Ak,)+F{k-Ak,)] . (42) 

Such deformations ensure that the mean value of 
the wave number k = k in the whole system and 
also the number of periodic units remains un- 
changed. 

For small values of Akj^ the expression at the 
right hand side of Eq. (|42] ) can be expanded in 
terms of a Taylor series and one obtains at lead- 
ing order in Akx'. 

F{k)=F{k)-^^-^^{Akxf + ... . (43) 

It depends therefore on the sign of the second 
derivative d^F{k)/dk^, whether the slight defor- 
mation of the periodic solution leads to a reduction 
or an enhancement of the free energy F{k) with re- 
spect to F{k). In the case of d^F{k)/dk^ < a si- 
multaneous small dilation and compression of the 
periodic solution of wave number k in neighbor- 
ing ranges leads to a reduction of the free energy 
with F{k) < F{k). I.e. for such parameter com- 
binations (£,^) periodic solutions are unstable. In 
the opposite case with d^F{k)/dk^ > a periodic 
solution with wave number k is stable with respect 
to longitudinal perturbations. Therefore, the curve 
separating in the e — k plane the range of stable 
from unstable solutions is determined via the con- 
dition 

3.4.2 Zig-zag stability boundary 

The zig-zag instability of a periodic solution of 
wave vector k is induced by perturbations with a 
wave vector perpendicular to k. Utilizing that the 
free energy of the perturbed solution with the wave 
vector k = (fc, Aky) only depends on the value of 
the wave number |k| = ^/k'^-\- (Aky)^, one finds 
for small values of Aky at leading order 

F(|k|)=F(fc) + l^(Afc/-f.... (45) 

In the range of k with dF{k)/dk < small 
transversal perturbations of the periodic solution 
reduce the free energy of the system and thus the 
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periodic solutions are unstable. A periodic solu- 
tion of wave number k with dF{k) /dk> is stable 
with respect to transversal perturbations. There- 
fore the zig-zag stability boundary in the £ — k 
plane is determined by the condition 



dF{e,k) 
dk 







(46) 



which corresponds for periodic solutions also to 
the minimum condition of free energy functional 
with respect to k. The conditions for the Eck- 
haus boundary given by Eq. (|44|) and the zig-zag 
line by Eq. (I46l) are valid for any two-dimensional 
isotropic system with a dynamics governed by a 
functional. 

Recipe for a determination of the stability 
boundaries of stationary periodic patterns in sys- 
tems with a potential: In a first step the nonlinear 
periodic solution of wave number k is determined 
either analytically (e.g. by a one-mode approxi- 
mation) or numerically. In a second step the free 
energy functional F{e,k) is determined for the 
periodic solution and in a third step the location 
of the Eckhaus stability boundary, £E{k), is deter- 
mined via the condition given by Eq. (I44l) and the 
zig-zag stability boundary by Eq. (|46] ). 

3.5 Stability in the weak segregation 
regime 

With the one-mode approximation given by 
Eq. (|40|) the functional in Eq. (fT4l) can be easily 
evaluated and the resulting free energy per period 
A =2ii/k\s given by: 



- Pb{k) _ 3 4 



(47) 



In this case the condition (l44l) leads to the follow- 
ing Eckhaus stability boundary: 



eE{k) 



3k^ + 5a^ 
F(F+3a) 



(48) 



This formula reads in terms of the reduced con- 
trol parameter r and the reduced wave number k as 
follows: 



rE 



3P-2^6-6P + 5 



(49) 



The zig-zag instability condition (|46l ) provides the 
wave number kzz (resp. kzz) at the zig-zag stabil- 
ity boundary: 



a 



1/4 



for k 



zz 



1] 



(50) 



This wave number does not depend on the control 
parameter e, respectively r. The lamella period 
Xe = 2% I kzz at the minimum of the free energy 
and the corresponding free energy per period #e = 
^(^zz) are given by 



1% 

"174 



and 



2v^) 



a 



(51) 



These results are obtained for a one-mode ap- 
proximation of the nonlinear periodic solution and 
an analysis of the related free energy agrees with 
a conventional stability analysis in terms of am- 
plitude equations, as described for the present 
system in more detail in Appendix |Bj 

The stability properties of periodic solutions 
of Eq. (fT2l) beyond threshold are summarized in 
Fig. [U The dot-dashed line represents the neu- 
tral curve r^iji) as described by Eq. (|28T ). In the 
range beyond r^(^ periodic solutions are in one 
spatial dimension only stable within the Eckhaus 
stability boundary, which is marked in Fig. [I]by 
triangles (full numerical analysis). For compari- 
son the Eckhaus-stability boundary in terms of a 
one-mode approximation is given by the dashed 
line, cf. Eq. (|49l ). The zig-zag stability bound- 
ary obtained by full numerical analysis is marked 
by the open circles in Fig. [Hand for the one-mode 
approximation in Eq. (l50l) by the solid line. To 
the left of the zig-zag boundary and to the right 
of the right Eckhaus boundary, spatially periodic 
solutions are unstable in two spatial dimensions. 

In a full numerical analysis the nonlinear pe- 
riodic solution is determined numerically, as de- 
scribed in more detail in Appendix |A] Then the 
stability boundaries are determined by the condi- 
tions given in Eq. (I44l) and in Eq. (|46l) , similar as 
for the one mode approximation above. 

Alternatively, the stability of nonlinear periodic 
solutions is determined by a linear stability analy- 
sis, as described in Appendix |A] and in Appendix 
|B] analytically in the range r > 0. The analysis of 
the functional and the linear stability analysis give 
identical stability boundaries. 
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Figure 1: The stability diagram of the periodic 
solutions of Eq. (fT2l) shows the neutral curve (N, 
dot-dashed line) according to Eq. (|28l) . the Eck- 
haus stability boundary obtained for the one-mode 
solution via Eq. (|49l ) (E, dashed line) and numer- 
ically (triangles), as well as the zig-zag stability 
boundary for the one-mode approximation given 
by Eq. (|50l ) (ZZ, solid line) and as obtained numer- 
ically (circles). The wave number km = y/l+r at 
the maximum of the growth rate o is given by the 
dotted line. 



The reduced control parameter range, < r < 1 , 
corresponds to a moderately deep quench of the di- 
block copolymer melt and belongs to the so-called 
weak segregation regime. In this range the Eck- 
haus and zig-zag boundary, as obtained by the one- 
mode approximation in Eq. (I40l ). coincide well 
with the numerical stability analysis, where sev- 
eral modes of the Fourier expansion in Eq. <B9\ 
are taken into account. In the range r > 1 the spa- 
tial shape of the nonlinear periodic solutions be- 
comes increasingly anharmonic and deviates from 
the one-mode solution given by Eq. (|40l ). Accord- 
ingly, the results of the one-mode approximation, 
as given by Eq. (|49l ) and Eq. (|50l) . start to devi- 
ate from the full numerical results obtained for the 
Eckhaus boundary (triangles) and the zig-zag line 
(circles). 

The dotted line in Fig. [T] shows the control pa- 
rameter dependence of the wave number k,n of 
the fastest growing perturbation with respect to 
the homogeneous state. The curve k,n{r) crosses 
the Eckhaus boundary at ^ = 5^/"^ ~ 1.5, corre- 



sponding to r = v^— 1 ~ 1.24. Therefore, after 
a deep quench with r > — 1 , the wave num- 
ber of lamellar structures developing during the 
early stages of microphase separation may lie in 
the unstable range to the right of the Eckhaus sta- 
bility boundary. The processes required to relax 
the wave number of the periodic solution back to 
the stable wave number band leads to the appear- 
ance of many defects as discussed further in Sec. |5] 

3.6 Strong segregation regime 

With increasing of the control parameter r (resp. 
e) the spatially periodic solutions of Eq. (fT2l) be- 
come rather anharmonic and for large values of r 
they can be approximated by a square wave of the 
form 

Y{x) = ao [ 1 - 20(jc - a /4) + 20(x - 3A /4)] , (52) 

with < X < A and the Heaviside step function 
&{x). In this limit the Green's function in Eq. (fT4l) 
is given by G{x,x') = \x — x'\/2. The gradient 
square term in the free energy Eq. ([HJ can be 
calculated by using the hyperbolic tangent profile 
x^f = =i=aotanh[(x — Xj)/^] with a small width ^ of 
the interfaces at xi = A/4 and X2 = 3A/4 as an 
approximation of the step function. Within these 
approximations one obtains the following expres- 
sion for the free energy per period A : 



j>j = — ^ £- — 



an 16 aX^ 



2 3A^ 48 



(53) 



A minimization of this expression with respect to 
ao gives 



al = £ 



16 aX^ 



3A^ 48 



(54) 



and with this amplitude the corresponding free en- 
ergy can be further simplified to 



(55) 



The equilibrium period of lamellae Ag correspond- 
ing to the minimum of the free energy is found by 
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a minimization of Eq. ( |55] ) with respect to A: 



1/3 



(56) 



In the strong segregation regime, i.e. Ec (large 
r), or equivalently for a —t- at a fixed e, the width 
of the interface ^ a/2/£ becomes small. In 
this case Xe at the minimum of ^ and the corre- 
sponding free energy per period are given by 



4(28] 
1 

~4 



1/6^-1/3 



(2ea; 



,1/3 



(57) 



The scaling Xe ~ cc^^^^ has to be compared with 
the scaling in the weak segregation limit Xg ~ 
a^i/"^ according to Eq. (ISTT ). The results for Ae 
and obtained in weak and strong segregation 
regime, which are given by Eq. (ISTT ) and Eq. (|57| |. 
respectively, are in each regime in good agreement 
with the results according to the full numerical so- 
lutions, as can be seen in Fig. [2l The transition 
from the weak segregation regime in the range 
a > 0.02 to the strong segregation regime in the 
range a < 0.02 is clearly visible in Fig. |2l The 
scaling exponents agree with those obtained ear- 
lier from variational calculations-^^ and numerical 
simulations .-^2^ 

4 Orientation of lamellae be- 
tween substrates 

The free energy of block copolymer films between 
two confining substrates depends on the orienta- 
tion of lamellae with respect to the substrates, on 
the distance between the confining substrates and 
on the surface properties of the bounding sub- 
strates. In order to model for instance a BCP 
film with the lamellae perpendicular to a substrate, 
which is in addition laterally confined by parallel 
side walls as in Refs.,— >^ we consider in this sec- 
tion Eq. (fT2l) in two spatial dimensions in the x — y 
plane with boundaries at y = and y = Ly. This 
analysis applies as well to BCP confined between 
two extended plane parallel substrates. 
The film thickness Ly is given in terms of the 



A. 
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Figure 2: The lamella period Xg at the minimum 
of ^{k) is shown in (a) as a function of a for e = 
1 and in (b) the corresponding free energy ^g ~ 
^{Xg). The solid lines are obtained by Eq. (ISTT ) 
(weak segregation), the dashed lines by Eq. (|57| ) 
(strong segregation). The bullets are the results of 
the full numerical calculations. 



dimensionless number d and the wavelength Xg'. 

Ly = dXg. (58) 

Boundary conditions. At plane substrates the 
boundary conditions for Xjr given by Eq. (fTSi) and 
Eq. (I2TI) take the following form. 



d\l/ 



y=0 



y=Ly 



0, (59a) 
(59b) 

g{WL,-W)\y=Ly, (59c) 



y=0,Ly 
giV-V0)\y=0 
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with Xj/o = Xj/siy = 0) and Xj/t, = Wsiy = Ly). 

Substrates preferentially wetted by one block of 
an AB copolymer are described by finite values of 
i//o and corresponding to so-called selective 
boundary conditions. We consider either symmet- 
ric selective boundary conditions at the two con- 
fining substrates, i//o = V^l,. 7^ 0, or antisymmetric 
ones, i//o = —WLy 7^ 0. Substrates being equally 
wetted by the A- and the B block of a copoly- 
mer correspond to neutral boundary conditions, 
V^o = Wl,. =0. As a third example of confined 
copolymer films we investigate also mixed bound- 
ary conditions, when one substrate acts like a se- 
lective boundary and the opposite one like a neu- 
tral boundary. 

With the wave vector = [k, 0) we describe the 
periodic order of lamellae perpendicular and with 
k|| = {0,k) the periodic order of lamellae parallel 
to the substrates. 

Numerical method: To find stationary solutions 
of Eq. (fT2l) with the boundary conditions Eq. (|59] ) 
a central difference approximation of the spatial 
derivatives is used. In the case of an orientation of 
lamellae parallel to the substrates one has to con- 
sider only the y dependence of Eq. (fT2l) and New- 
ton's iteration method is used for its solution. For 
lamellae perpendicularly oriented to the substrates 
two-dimensional simulations of Eq. (fT2l) are re- 
quired. In this case we use a simple relaxation 
method with the width of one period = along 
the X direction. 

For a given solution ^f{r) the total free energy 
F = F}, -\- Fs h calculated by integrating Eq. ([14] ) 
and Eq. (fT9l) numerically. In order to determine 
the last term in Eq. ([T4l) . Poisson's equation for 
the auxiliary function in Eq. (fT6] ) is solved numer- 
ically by a relaxation method. The spatial dis- 
cretization was chosen to be dx = 5y = 0.5 for 
most of the calculations, which provide a relative 
error of the free energy less than 1%. For a transi- 
tion range, where the free energies of lamellae par- 
allel and perpendicular to the substrates become 
comparable, the discretization was decreased to 
5x = 5y = 0.25 for the purpose of a precision 
higher than 0.2%. Note, that the values for the 
total free energy F presented in the following are 
divided by the system size A = L^x Ly, i.e. we use 
the free energy per unit system size. 



4.1 Selective boundary conditions 

In the case of homogeneous, selective bound- 
ary conditions with i//o, Vl, 7^ at the substrates, 
lamellae parallel to the boundaries have a lower 
free energy than perpendicularly oriented ones, as 
shown in this section. If the values i//o, V/l,, have 
a magnitude similar to the maximum of the ampli- 
tude of \]/{y) in the bulk, then the envelope of \f/{y) 
is only slightly deformed near the boundaries. 
In the case i//o and 1//^,, agree with the extrema of 
in the bulk, then also the boundary condition 
d^'Wl n, ~ can be fulfilled at an extremum of 

> ^ \y=0,Ly 

a periodic function without any deformation. 




y/K. y/K 



Figure 3: Four stationary solutions of 
Eq. (fT2l) in a film of thickness Ly = 5A,e are shown 
for symmetric selective boundary conditions i//o = 
V/l, = 1- Part (a) and (c) show periodic solutions 
with an integer number of lamellae and part (b) 
and (d) so-called (unstable) saddle-point solutions. 
Parameters e = g = I and a = 0.015. 

Examples of \]/{y) for lamellae parallel to the 
boundary are shown in the case of selective bound- 
ary conditions, i//o = V^l, = 1, in Fig. ^a) and 
Fig. [3]^c) for a copolymer film of thickness Ly = 
5Xe in the strong segregation regime. The periodic 
field in Fig. ^a) and in Fig. [3tc) differs in 
the number of periods on the interval [0,Z^], cor- 
responding to different values of the wave number 
k. The wavelength of the solution with five periods 
in Fig.[3];c) corresponds to X = Xe at the minimum 
of the free energy. In Fig.|3]^a) the solution has four 
periods with a wavelength X > Xe and this station- 
ary solution in an unconfined system is unstable 
according to the results in Fig. [1] However, dy- 
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namical simulations show, that this wavelength is 
stabilized in a confined thin film. The solution in 
Fig. life) has six periods on the interval Ly = SXg 
and a wavelength A smaller than Xe. This solu- 
tion is according to the results presented in Fig. [1] 
expected to be stable in unconfined systems. The 
free energy F /Ly of the solutions in Fig. Ufa) and 
(e) is larger than in Fig. life) at the minimum of the 
free energy. 
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Figure 4: The normalized free energy F /Ly oi 
a copolymer film of thickness Ly = SAg is shown 
as a function of the normalized width A / (solid 
line) for lamellae parallel to selective boundaries 
with i//o = M^Ly = 1 and parameters as in Fig. |3] 
(a)-(d) mark the free energy of the corresponding 
solutions in Fig. |3l 

The stationary solutions in Fig. |3]^b) and 
Fig. ^d) are so-called saddle point solutions, 
which are unstable. As the characteristic wave- 
length A of these two solutions we take the dis- 
tance between two extrema in the "undistorted" 
range of each solution. With this definition the 
solution in Fig.[3td) has a wavelength between the 
wavelengths of the two solutions in Fig. Oe) and 
(c) and the saddle point solution in Fig. [3tb) has a 
wavelength between that of the periodic solutions 
in Fig.[3ta) and (c). The free energy of both saddle 
point solutions is higher than that of the periodic 
solutions marked as (a), (c) and (e) in Fig. |4l The 
locally strong deformation of the periodic solu- 
tions in Fig. [3tb) and (d) may occur at different 
locations y in the region (0,Ly), depending on the 
initial profile. In order to include or to remove one 
periodicity, as for instance by changing from the 
solution given in Fig. [3]^ a) to Fig.[3]^c) or reversely, 
the local maximum of the saddle point solution 



Fig. [Stb) has to be crossed. Such energy barriers 
are essentially responsible that states with a wave 
number k ^ kc are stable in BCP films even if the 
wave number does not correspond to the minimum 
of the free energy. 

1 I ' ' ' ' 1 



0.5 




y/K 

Figure 5: Stationary solutions of Eq. (fT2] ) for 
lamellae oriented perpendicularly to the bound- 
aries at a position xi (dashed), where ^f{x,y) takes 
its maximum in the bulk and at xi (solid), where 
^f{x,y) takes its minimum. The same parameters 
as in Fig.[3l 

The y dependence of the order parameter is 
rather different in the case of the lamellae perpen- 
dicular to selective boundaries. In this case ^f{x,y) 
is a periodic function along the x direction and 
selective boundary conditions force a finite value 
V^o,Lv 7^ (Vb.Lv > 0) at ); = 0,Ly, independent of 
the phase of the function. At positions xi, where 
^f{x,y) takes its maximum in the bulk, the order 
parameter i//(xi , y) is nearly undeformed as a func- 
tion of y, as can be seen by the dashed line in 
Fig. [51 However, the imposed selective boundary 
condition requires strong deformations of i//(jC2,j) 
along the y direction at positions X2, where ^f{x,y) 
takes its minimum in the bulk, as indicated by the 
solid line in Fig. [51 As a consequence of such 
strong deformations of i//(jC2,y) near the bound- 
aries, perpendicularly oriented lamellae have for 
selective boundary conditions a higher free energy 
than parallel oriented ones as shown in more detail 
in Fig. [61 

The free energy per unit area, Fy/A, of a copoly- 
mer film with its lamellae parallel to the bound- 
aries has as a function of the film thickness Ly ~ 
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Figure 6: The free energy per unit area F/A as. 
a function of d for lamellae parallel or perpendic- 
ular to symmetric selective boundaries in (a) with 
V^o = Vty = 1 and in (b) with i/Aq = V^l, =0.5. Pa- 
rameters a = 0.015 and e = g = 1. 

dXe local minima at integer values of d as indi- 
cated by the solid lines in Fig. [6ta) and Fig. MX>)- 
For parameters used in Fig. |6la) the corresponding 
minima of the solid line have even an equal height. 
The dashed line in Fig.[6ta) shows the normalized 
free energy, F±/A,of lamellae perpendicular to the 
substrate, which is for nearly all values of d higher 
than F^^/A. The periodically occurring strong vari- 
ation of the order parameter i//(x2,3') for perpen- 
dicularly oriented lamella in the case of selective 
boundary conditions, as indicated at one position 
X2 by the solid line in Fig. [51 enhances the free en- 
ergy compared to the nearly undeformed function 
^^{y) in Fig. 12 a) for parallel lamellae. 

The decay of F±{d) /A in Fig.[6ta) indicates that 
the weight of the strong deformation of \l^{x,y) 
of perpendicularly oriented lamellae near the sub- 



strate becomes smaller with increasing thickness 
of the copolymer film. In Fig. [6ta) only for a very 
thin film-thickness of about <i ~ 1 .5 the free energy 
of parallel oriented lamellae is higher than for per- 
pendicularly oriented ones. 

In Fig. [6tb) the normalized free energies of 
parallel and perpendicularly oriented lamella are 
shown in the case of a reduced selectivity, i//o = 
V/l, = 0.5. Since the control parameter e (resp. 
r) is unchanged compared to the case in Fig. |6ta), 
a reduced value i/Zq = V^l, =0.5 requires now a 
deformation of the function at the bound- 

aries also for parallel lamellae. This deforma- 
tion increases the normalized free energy of par- 
allel oriented lamellae, while the normalized free 
energy of perpendicularly oriented ones remains 
nearly unchanged, as can be seen by comparing 
the dashed lines in Fig.[6ta) and (b). This enhance- 
ment of the free energy is stronger for small values 
of d than for larger values of d, because of the de- 
creasing weight of boundary effects with increas- 
ing film thicknesses. 

As a consequence of this energy enhancement in 
the case of a reduced preferential adsorption, there 
are now two maxima of the free energy of parallel 
lamellae in Fig.[6tb), at about d^ 1.5 and J ^ 2.5, 
where the free energy is higher than that of lamel- 
lae perpendicular to the substrates. Such situa- 
tions of confined diblock copolymers were also 
studied experimentally in thin films in the range 
d= 1 .4 — 3.2 by varying the selectivity of the sub- 
strates. Here, a reduction of the preferential ad- 
sorption leads at about d ~ 2.5 to a frustration and 
lamellae perpendicularly oriented to the bound- 
aries. However, in agreement with our simula- 
tions, for block copolymer films being less frus- 
trated and also for strong preferential adsorption 
the parallel orientation of lamellae remains always 
preferred in this experiment. 

The normalized free energy of lamellae parallel 
to selective substrates becomes smaller with de- 
creasing values of a, as shown in Fig.|71 This trend 
is similar to the a-dependence of the bulk free en- 
ergy given by Eq. (|5TI) [see Fig. l2lb)]. Since de- 
creasing values of a correspond to increasing val- 
ues of the thickness of the lamellae, the weight of 
surface effects decreases, that leads to a reduction 
of the peak height with a, as can be seen in Fig. |7] 
too. 
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Figure 7: The free energy per unit area F/A of 
lamellae parallel to the substrate as a function of d 
for selective boundaries \]/q = x^^. = 1 and differ- 
ent values of a and £ = s= 1 . 



For asymmetric selective boundary conditions at 
the substrates, when one of the two substrates is 
preferentially wetted by one block and the other 
one by the second block of the copolymer, the nor- 
malized free energy of parallel lamellae has local 
minima at a film thickness close to a half-integer 
multiple of the equilibrium lamellar thickness A^, 
i.e. for d = 1.5,2.5,3.5,..., as indicated by the 
solid line in Fig. [H A situation with compara- 
ble free energies for lamella orientations paral- 
lel and perpendicular to asymmetric boundaries is 
only met in the range of very thin films of about 
d ^ 2. Otherwise the trend, that lamellae paral- 
lel to the substrates have for asymmetric selective 
boundaries a lower free energy, can be explained 
by the same arguments as given above for the case 
of symmetric selective boundary conditions. 

By a reduction of the surface interaction strength 
g (leading to non-interacting or quasi-periodic 
boundary conditions for ^ — 0) or by a reduction 
of the preferred difference \j/s between the con- 
centrations of A- and B blocks at the boundary 
the free energies of both orientations can become 
comparable in the range of very thin films like in 
Fig. Ob). However, in the range of thick films a 
parallel orientation of lamellae is always preferred 
in the case of selective boundary conditions. De- 
tailed studies on that issue can be found elsewhere 
(see, e.g.,—"—). 
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Figure 8: The free energy per unit area F/A of 
parallel and perpendicularly oriented lamellae as a 
function of d for asymmetric selective boundaries: 
V^o = 1, V^Ly = -1- Parameters a = 0.015 and e = 
5=1- 

4.2 Neutral boundary conditions 

Neutral boundaries with Yq = Xj/^^, = correspond 
to substrates, which are neither preferentially wet- 
ted by the A- nor by the B block of a copolymer. 

A first estimate of the expected preferred lamel- 
lae orientation may be gained by considering the 
effect of neutral boundaries in the weak segrega- 
tion limit with small values of r > 0. In this range 
a representation of i// = A exp (? k|| ^ ■ r) + cc as in 
Eq. (I32l) is useful, where = (A:, 0) is the wave 
vector in the {x,y) plane of lamellae perpendicu- 
lar and k|| = (0,fc) of lamellae parallel to the sub- 
strates. The envelope A{x,y) decays in the case 
of neutral boundaries from its finite bulk value 
A oc yfp to the boundary value A ^ 0. Such a re- 
duction of the envelope causes an enhancement of 
the free energy per unit size compared to the case 
without boundary effects. 

The transition layer, in which the envelope 
A{x^y) changes from its bulk value to that at the 
boundary, is for perpendicularly oriented lamellae 
according to Eq. (|38T ) proportional to ^2 °^ r^^^^, 
which is for small values of r smaller than the 
transition layer of parallel oriented lamellae pro- 
portional to (^1 oc r^^/^. Since the transition range 
is smaller in the case of lamellae perpendicular 
to the boundaries, we expect a smaller energy of 
perpendicularly oriented lamellae than for parallel 
oriented ones. 
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Figure 9: Stationary solutions of Eq. (fT2] ) for neu- 
tral boundary conditions VA) = Wty = 0: In (a) at 
a position x\ where takes its maximum 

for lamellae perpendicular and in (b) parallel to 
the boundaries. In (b) the dots mark the analyti- 
cal approximation as described in Appendix O for 
the same boundary condition. The parameters are 
Ly = lOAe, g = \,r = 0.021, k = 0.7 (correspond- 
ing to £ = 1 and a = 0.24). 



Full numerical solutions of Eq. (fT2l) by taking 
into account the boundary conditions (|59| ) with 
Vb = Vl, = are shown in Fig.|9]in the weak seg- 
regation limit at r = 0.021 for perpendicularly ori- 
ented lamellae in part (a) and for parallel oriented 
lamellae in part (b). In Fig.|9tb) we show also the 
analytical approximation of the solution (symbols) 
for the same boundary conditions, as described in 
the Appendix O 

As indicated by the estimate in the previous 
paragraph, the length of the envelope of ^/(j) 
needed for the transition from its value in the bulk 
to that at the boundary is indeed larger for par- 
allel oriented lamellae than for the perpendicu- 



larly oriented ones. A narrower transition range 
causes a smaller enhancement of the free energy 
and therefore, in the range of small values of r > 
(weak segregation limit) lamellae perpendicularly 
oriented with respect to the substrates are ener- 
getically preferred. This behavior extends also to 
the strong segregation regime with larger values of 
r, as we have tested by further numerical calcula- 
tions. 
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Figure 10: The free energy per unit area F /A of 
parallel and perpendicularly oriented lamellae as a 
function of d for neutral boundaries i//o = M^Ly = 
and the parameters a = 0.015, £ = ^ = 1. 

For numerical stationary solutions of Eq. (fTZl i in 
the strong segregation regime the free energy of 
lamellae, that are perpendicularly oriented to neu- 
tral boundaries, does not differ very much from 
the free energy obtained in the case of selec- 
tive boundaries, as can be seen by comparing the 
dashed curves in Fig. |6] and Fig. \T0\ On the other 
hand, the decay of the envelope of \ff{x,y) close 
to the boundaries, as shown in Fig. Hfb), enhances 
the free energy of parallel lamellae compared to 
the case of selective boundaries, cf. Fig. |6l As a 
consequence of both trends, in the case of neutral 
boundaries perpendicularly oriented lamellae have 
always a lower free energy than parallel oriented 
ones, as shown in Fig. \T0\ 

Note, that the free energy of parallel oriented 
lamellae has in the case of neutral boundaries local 
minima as a function of the film thickness close 
to integer and half-integer multiples of A^,. The 
tendency of lamellae to align perpendicularly to 
the substrates in the case of neutral boundaries has 
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been also found in Refs.— 

In the weak segregation limit, i.e. small values of 
r, approximate analytical solutions of Eq. (fT2l) are 
derived for lamellae parallel to the substrates, as 
described in more detail in Appendix O Depend- 
ing on 1^0 and xif^, such an analytical approxima- 
tion can be very good as can be seen for example 
in Fig.Hb). 

4.3 Selective versus neutral bound- 
aries 

It depends on the ratio between the extremal val- 
ues of the amplitude of i// in the bulk and the in- 
duced values at the boundaries whether the bound- 
ary conditions act more like selective or neutral 
boundary conditions. This can be recognized for 
instance by comparing the difference between the 
free energy of parallel and perpendicularly ori- 
ented lamellae in Fig. [6] and Fig. [10] for the three 
different values: Xj/Q = Xj/^, = 1, Xj/o = Vz-y = 0.5 
and y/Q = yr^, = 0. While in the case y/Q — y/^, = 1 
the maximum of yf{y) in the bulk is similar to the 
imposed value at the boundary, in the other two 
cases the maximum in the bulk is larger than at the 
boundaries. 

The ratio between the maximum value and the 
value at the boundary can also be changed by 
changing the quench depth, i.e. by changing r (re- 
spectively £), but keeping now the values i//o = V^l,, 
fixed. In this case the maximum bulk value can be 
either smaller or larger than the values at bound- 
aries, depending on r. We found that for different 
values of i/Zq = V^l, variations of the control param- 
eter r do not induce a reorientation of the lamellae 
with respect to the boundaries. 

4.4 Mixed boundary conditions 

The results presented in the previous sections in- 
dicate that a combination of a selective and a neu- 
tral boundary condition may lead to almost equal 
energies for parallel and perpendicularly oriented 
lamellae over a large range of values of the film 
thickness Ly,. Therefore, we compare in this sec- 
tion for mixed boundary conditions the free ener- 
gies of homogeneously oriented lamellae parallel 
and perpendicular to the confining substrates. 
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Figure 11: The free energy per unit area F/A of 
parallel and perpendicularly oriented lamellae as a 
function of d for mixed boundary conditions: In 
(a) \j/o = 0, \j/L, = 1 and in (b) y/Q = 0, y/t, = 0.5. 
Parameters a = 0.015 and £ = ^ = 1. 

In Fig. [TT] the free energy per unit size is shown 
as a function of d for perpendicularly (dashed 
lines) and parallel (solid lines) oriented lamellae in 
the case of mixed boundaries with i/Zq = and ei- 
ther yr^^, = 1 or yr^. =0.5. One may compare these 
results with those given in Fig. |6] for symmetric 
selective boundary conditions yro = yf^, = 1, 0.5. 
There are two major differences between the re- 
sults in both figures. The energy differences be- 
tween the two lamellae orientations are smaller for 
mixed boundary conditions in Fig. [TT]and the free 
energy of parallel oriented lamellae now has local 
minima at integer and half integer values of d. 

The trend indicated in Fig. \TT\ suggests that by 
changing the boundary condition at one surface 
from selective to neutral and keeping the other sur- 
face neutral, the preferred lamellae orientation can 
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be changed from parallel to perpendicular. This is 
0.02 




0.4 0.6 



Figure 12: The relative free energy difference 
{F^ — F[| ) / I is shown for three thicknesses Ly = 
3Ag, 6Ae, lOXg of the BCP film as a function of 
the selectivity Xj/^^. at one boundary and i//o = at 
the opposite one. Beyond a critical value \if^{crit) 
the preferred lamella orientation is parallel to the 
boundaries and perpendicular below. Parameters 
a = 0.015 and e = g=l. 

shown in Fig. [I2l where the relative difference of 
the free energy {F± —Ft^y)/\F± \ is plotted as a func- 
tion of 1//^, for three different values of the film 
thickness Ly = dXe- Fig. [T2l shows in addition that 
the critical value \j/^,{crit), where both lamellae 
orientations have the same free energy, is rather 
independent of the film thickness d. This may be 
explained as follows. For the parameters used in 
Fig.[T2lthe two length scales introduced in Eq. (|38T ) 
are nearly equal and both are small: i^i/Ag ~ 0. 13 
and ^2/ K ^O.l. I.e. the influence of the boundary 
is similar for the three values of the film thickness 
in Fig. [T2I only the weight of the influence is re- 
duced by increasing the film thickness. The prior 
effect leads to a smaller slope of the curves with 
larger values of d in Fig.[T2l 

The critical selectivity \if^{crit) depends weakly 
on the parameter a. For smaller values of a and 
therefore a larger lamella period Ag, the critical 
value of \if^.{crit) is smaller. Thus lamellae with 
a larger period require a smaller selectivity of the 
surface to realign. 

Note that in case of relatively thin films the free 
energy of parallel oriented lamellae as a function 
of thickness shows pronounced oscillations be- 



tween local minima and maxima [see Fig. [TTT a)1. 
This leads to a weak thickness dependence of 
the critical selectivity when considering the thick- 
nesses that correspond to a maximum and a mini- 
mum of i^i (see Fig.[T3l). In the case of a maximum 
of i^i the reorientation takes place at higher values 
of x^Ti^, than for a minimum. With increasing film 
thicknesses this difference is rapidly decreasing. 




Figure 13: The relative free energy difference 
(F^ — F^y)/\F^\ is shown as a function of the se- 
lectivity at one boundary with \]/q = at the 
opposite one for two BCP films with Ly — IJSXg 
(lower curve) and Ly = 3Ae (upper curve). In both 
cases the preferred lamella orientation is parallel 
to the boundaries beyond a critical value \if^,{crit) 
and perpendicular below. Parameters a = 0.015 
and e — e—l. 



5 Dynamics of microphase sep- 
aration 

The spatio-temporal dynamics of microphase sep- 
aration in copolymers in two spatial dimensions 
between two parallel boundaries and the related 
lamellar (orientational) order is investigated here. 
We also describe typical differences between the 
evolution of structures in the strong and the weak 
segregation regime in BCP films confined between 
different boundaries on the one hand and in uncon- 
fined systems on the other hand. 

For numerical simulations of Eq. (fT2]) we use 
a central difference approximation of the spatial 
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derivatives with dx = dy = 0.5 and an Euler inte- 
gration of the resulting ordinary differential equa- 
tions with a time step At = 10^^ — 10^-^. In the un- 
confined case periodic boundary conditions are ap- 
plied and a system size = Ly = 256 (~ 14A<j for 
a = 0.015) is chosen. For block copolymer films 
of thickness Ly = 6Xe between two substrates, dif- 
ferent combinations of the boundary conditions 
along the y direction are used, cf. Eqs. (|59] ). and 
periodic boundary conditions along the x direc- 
tion with = 4Ae, 8Ae, or 32Ae. To mimic a 
quench we start simulations of Eq. (fT2] ) with ran- 
dom initial conditions for i// of a small amplitude 
of about 10^^. Typical scenarios of the dynamics 
of microphase separation are studied in the strong 
segregation regime at a control parameter £ = 1 
(r = 3.08) and in the weak segregation regime at 
£ = 0.37 (r = 0.5). 

Microphase separation can be characterized by 
the structure factor of the evolving patterns. 



with 



s{Kt) 



e'^-'\l/{r,t)dr . (60) 



Since we expect anisotropy effects in BCP films 
confined between two substrates we introduce dif- 
ferent characteristic lengths along the x and the y 
direction: 

l,it)=7t/{k,){t), ly{t) = n/{ky){t). (61) 

The averaged wave numbers (fc,) are 

r^max rAk 

/ dkx / dky Si^kxtkyft^kx 
= ■ , (62a) 

10 J-Ak ' 




rAk 

dky I dk^ S (^kx 7 ky ^ ^ ) ky 

m) = " ^'.M ' ' (62b) 



dky dkx S{kx,ky,t) 

" J-Ak 



where Ak is the half-width of the corresponding 
peak of the structure factor along kx and ky, re- 
spectively. 



5.1 Unconfined systems 

During microphase separation in diblock copoly- 
mers the most unstable perturbation with respect 
to the homogeneous basic state i// = has the wave 
number k,„ — y^£/2, cf. Eq. (|271) . similar as in bi- 
nary mixtures.— The coarsening regime in diblock 
copolymers below Tc is at early stages similar as in 
polymer blends, as indicated by two snapshots of 
a simulation in Fig. [141 In diblock copolymers the 
coarsening process of phase separation is limited 
by the chemical bond between an A- and B block, 
which limits the domain size of phase separation 
to the order of the chain length of diblock copoly- 
mers. 




(c) 



(d) 



Figure 14: Microphase separation in the strong 
segregation regime for the parameters £ = 1 
(r = 3.08) and a = 0.015 is shown in a two- 
dimensional system with Lx = Ly. = 256 ~ 14Ae 
(Xe = 19.3) and periodic boundary conditions: In 
(a) at the time t = 10^ with the average wave- 
length of about A ^ 11.0 and in (b) at t = 10"^ 
with A ~ 17.5. Dark and bright regions in the top 
part correspond to A- and B block rich phases, re- 
spectively. The bottom part (c) and (d) show the 
corresponding structure factors, where the bright 
regions correspond to large values of 5(k, t) . 



The pattern shown in Fig. \T4t h) has an average 



18 



wavelength still below the optimal wavelength 
at the minimum of the free energy. With a further 
progress of time the mean wavelength approaches 
only very slowly towards Xg, because the system 
has to get rid of lamellar imperfections by diffu- 
sion processes. 
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Figure 15: The temporal evolution of the charac- 
teristic length scales /x(?) and ly{t) (averaged over 
10 independent runs) after a quench is shown in 
the strong segregation regime for the same param- 
eters as in Fig. [14] and periodic boundary condi- 
tions. 

In an unconfined system the two length scales 
lx{t) and ly{t) evolve in a similar way and this 
isotropic behavior of the block copolymer melt is 
also reflected by the rotational symmetry of the 
structure factor shown by the parts (c) and (d) of 
Fig. \T4\ During the intermediate regime between 
the early stage of phase separation with a domi- 
nating growth of the perturbation of wave number 
k,„ and the late stage of coarsening with an aver- 
age domain size, ~ ly ~ Ae/2, one observes the 
scaling ~ ly ~ t^^^ as shown in Fig. \T5\ which 
is common for polymer blends. Such a scenario is 
typical for a deep quench into the strong segrega- 
tion regime at about £ ~ 1 and beyond. 

In the weak segregation regime (e = 0.37 or 
smaller) the wave number of the fastest growing 
mode during the early stage of phase separation, 
k,„ = \/£/2, is already closer tolTt/Xg and there- 
fore, one observes during pattern evolution only 
a slight variation of the scales ljf{t) and ly{t), as 
shown in Fig. [161 Typical snapshots of phase sep- 
aration during the early stage in the weak seg- 
regation regime are similar to patterns shown in 
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Figure 16: The temporal evolution of the charac- 
teristic length scales lx{t) and ly{t) (averaged over 
10 independent runs) after a quench is shown in 
the weak segregation regime at £ = 0.37 (r = 0.5) 
for a system with periodic boundary conditions 
and otherwise the same parameters as in Fig. [14] 
and in Fig. [T5l 

Fig.[l4l:b). 

5.2 Confined systems 

For a block copolymer film confined between two 
selective boundaries with i//o = V^l, = 1? three 
snapshots during microphase separation in the 
strong segregation regime at £ = 1 (r = 3.08) are 
shown in Fig. [ITl which will be compared later 
with results in the weak segregation limit. 

After a deep quench at ? = far below threshold 
£c = (for a = 0.0015 one has £,. = 0.077) 
the selective boundary conditions trigger close to 
the substrates immediately a large value of the or- 
der parameter i// and in this case lamellae orient 
parallel to the substrates. One should note, that 
for the parameters in Fig. [T7]the wavelength of the 
fastest growing mode A,,, = 2n/k,n = 2n^/2/e — 
2%^2 ^ 8.9 is much smaller than the wavelength 
Xe = 4(2£)i/^a"i/^ ^ 39 at the minimum of the 
free energy of a lamellar structure. 

In the strong segregation regime the correlation 
lengths in the x and y direction are rather small and 
therefore a surface induced orientational order of 
the lamellae occurs only within short ranges near 
the boundaries as indicated by parts (a) and (b) 
in Fig. [17] Further away from the substrates the 
orientation of lamellae is only weakly influenced 
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Figure 17: Microphase separation in a BCP film 
between two selective boundaries with i/Zq = V^l, = 
1 is shown in the strong segregation regime at £ = 
1 (r = 3.08) in (a) at time t = 20 after the quench, 
in (b) at t = 40, and in (c) at t = 10^ close to the 
final state. Parameters = 4Ae, Ly = 6Xe, (X = 
0.0015 and ^= 1. 

by boundary conditions and the lamellae are disor- 
dered, whereby this disorder resembles very much 
to that in the unconfined case, as shown in Fig. [14] 
and has also been observed in experiments on con- 
fined thin films.— >^ The average wavelength of 
the structure tends in the long time limit to Xg. 

Typical lamellar structures at the late stage of 
microphase separation in the strong segregation 
regime are shown in Fig. [18] for three types of 
boundary conditions. Neutral boundary conditions 
¥o = Vz-v = are used in Fig.[T8];a), symmetric se- 
lective boundaries i//o = V^l,. = 1 in Fig.[T8lfb') and 
in Fig. \W[c) mixed boundary conditions, i/Zq = 0, 
V/l,. = 1 • The simulations of Eq. ([T2]) were started 
with random initial conditions. 

In the case of symmetric selective boundaries in 
Fig. \W[h) lamellae parallel to the substrates have 
the lower free energy as shown for the defect free 
lamellar order in Fig. [6] For neutral boundaries in 
Fig. [T8]' a) an orthogonal lamellae orientation close 
to a boundary is favored, which is in agreement 
with the results shown in Fig. \T0\ In the case of 
mixed boundary conditions in Fig. [TST c) lamellae 
are oriented parallel close to the selective (upper) 
boundary and perpendicular to the neutral (lower) 
boundary. 

The pattern away from the boundaries in the 
bulk shows for neutral boundaries [Fig. [Hla)] a 
stronger disorder compared to the case of selective 
boundaries [Fig. [TSlb)!. In the strong segregation 




Figure 18: Microphase separation is shown in 
the strong segregation regime in a confined sys- 
tem at ? = 10^ after a quench: In (a) for neu- 
tral, i//o = V^Lv = 0> in (b) for symmetric selective, 
V^o = V^L, = 1 > and in (c) for mixed boundary con- 
ditions, i/Zo = 0, 1//^,, = 1 . Parameters = 32Xe, 
Ly = 6Xe, a = 0.015, e = 1 (r = 3.08) and g = 1. 

regime the coherence lengths are small for both 
cases, but in the case of selective boundary con- 
ditions larger values of i// close to the boundary 
are induced and this causes a more regular lamel- 
lae orientation in the bulk. 

The pattern Fig. [TST b) consists of seven periods 
parallel to the horizontal x-axis around the center 
of the image and six periods with A = Ag close to 
the right end, whereby both regions are connected 
by a pattem including defects. The wave numbers 
corresponding to seven and six lamellae between 
the boundaries lie both in the range k/kc > 1 of 
the stability diagram in Fig.[T] where a straight and 
defect free lamellar order is linearly stable with re- 
spect to small perturbations. 

In simulations started with random initial condi- 
tions, patterns with k,n > kc grow with the largest 
rate and therefore a lamellar order with small 
wavelengths is preferred during the early stage of 
microphase separation. Since one has in Fig. [TJa 
wide wave number range k> kc of stable straight 
lamellae, a relaxation of a pattem like in Fig. [TST b) 
to the homogeneous state with six lamellae, which 
has the lowest free energy, is a long lasting pro- 
cess. 

We showed in the previous section in Fig. \TT\ a) 
that in the case of mixed boundary conditions and 
parameters as in Fig. [TSl c) a defect-free order of 
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lamellae parallel to the substrates has a lower free 
energy than perpendicularly oriented ones. In sim- 
ulations of extended systems with mixed bound- 
aries, where patterns with defects may occur, nei- 
ther a parallel nor a perpendicular orientation of 
lamellae is preferred in the bulk. Moreover, the 
pattern in Fig. [TSlc) shows, that for mixed bound- 
ary conditions an orientational transition across 
the block-copolymer film can be expected, from 
parallel oriented lamellae at the selective (upper) 
boundary to a perpendicular lamellae orientation 
at the neutral (lower) boundary. The free en- 
ergy of the pattern in Fig. [TSt c) is higher than the 
free energy of parallel oriented lamellae as shown 
Fig. Eta) and lower than that of perpendicularly 
oriented ones. 
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shown in Fig.[T9ta). is rather similar to the uncon- 
fined case shown in Fig. [151 During the early stage 
of phase separation the dominating length scale 
is again that of the fastest growing mode, which 
is followed by the intermediate coarsening regime 
with Ix ~ ly °^ t^^^, before lx{t) and ly{t) terminate 
again at the typical length scale Xe/2 of a diblock 
copolymer. 

In the case of symmetric selective boundary con- 
ditions, v/q = V^L, = 1' the length scales and 
ly exhibit a different behavior during the initial 
stage of phase separation and especially the be- 
havior of lx{t) is changed significantly, as shown 
in Fig. [T9lb). A comparison of Fig. \T7\ a) and 
Fig. WTl h) reveals, that during the early stage of 
microphase separation compositional waves are 
induced by the selective boundaries and they prop- 
agate into the copolymer film. These induced 
composition waves near the boundaries have a 
quasi-infinite wavelength along the x direction, 
while the wavelength ly along the y direction be- 
haves similar as in Fig. [Uta) for neutral bound- 
aries. Far away from the selective boundaries one 
finds a random lamellae orientation and therefore 
Ix behaves in the bulk at an intermediate and late 
stage of microphase separation similar as in the 
case of neutral boundaries. 



(a)l 



(c)l 
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Figure 19: The time dependence of the char- 
acteristic length scales lx{t) and ly{t) (averaged 
over three runs) is shown in the strong segregation 
regime for the same parameters as in Fig. [18] (a) 
neutral boundaries, i/Zq = WLy = 0, and (b) selective 
boundaries, i/Zq = Wty = 1- 

The temporal evolution of the lengths lx{t) and 
ly{t) in the case of neutral boundary conditions, as 



Figure 20: Microphase separation is shown in a 
confined system in the weak segregation regime 
with £ = 0.37 (r = 0.5) at the time t = 500 after 
a quench and for three boundary conditions: (a) 
neutral, (b) symmetric selective and mixed in (c). 
Other parameters as in Fig. [18] 

For comparison, we show in Fig. [20] late stage 
patterns in the weak segregation regime at £ = 0.37 
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(r = 0.5) for the same confined systems as in 
Fig. [T8l These patterns show at a time t = 500 
already a similar order as in the strong segregation 
limit at ? = 10"^ (Fig.[l8]) in spite of the fact that the 
dynamics is slower for smaller values of the con- 
trol parameter r. However, by a reduction of the 
control parameter from r = 3.08 in the strong seg- 
regation limit to r = 0.5 one has an enhancement 
of the length scales from ^ 0.13 to ^\ ^ 0.32 
and from ^2 ~ 0.10 to ^2 ~ 0.16. These higher 
coherence lengths increase simultaneously the ac- 
tion range of the boundaries and both effects cause 
a higher lamellar order in a thin copolymer film 
within a shorter time (Fig. |20] ). 

Similar as in the strong segregation regime 
one finds in the case of selective boundaries 
[Fig. |20lb)1 again a higher lamellar order than in 
the case of neutral boundaries [Fig. |20l a)1. This 
is in agreement with the observation, that is for 
r = 0.5 roughly by a factor of 2 larger than ^2- This 
reasoning also confirms the results for the case of 
mixed boundaries [Fig. |20tc)1. where the action 
range of the neutral (lower) boundary is smaller 
than that of the selective (upper) boundary. 

The characteristic lengths l^{t) and ly{t) develop 
for r = 0.5 and neutral boundaries again very sim- 
ilar as in the unconfined case in Fig. \T6\ In the 
case of selective boundaries [Fig. ISOl b)! the two 
scales lx{t) and ly{t) show a similar behavior as in 
the strong segregation limit, only the saturation of 
ly{t) takes place already at r ~ 10^. 

In summary selective boundary conditions are 
more efficient for controlling the orientation of 
lamellae in copolymer films than neutral ones. A 
comparison of the results in Fig. [Island in Fig. |20] 
suggests in addition that a quench to a small value 
of £ > £c-, followed by a further enhancement of 
£ into the strong segregation regime, favors a co- 
herent order of the lamellae. In the case of mixed 
boundary conditions one obtains "mixed" lamel- 
lar structures as shown in Fig. iTSTc). which can be 
also interpreted as a coexistence of two different 
boundary induced lamellae orientations. 

As described in Sec. HI different numbers of par- 
allel oriented lamella between selective substrates 
can have at certain values of the distance Ly be- 
tween the boundaries equal free energies. For ex- 
ample at Ly = 4.4373 and for parameters as given 
in Fig. |2T| solutions with five and four lamellae par- 



allel to the substrates have the same free energy. 
Such a coexistence is shown in Fig. [2T] where the 
interface between both solutions does not move. 
This coexistence has a strong similarity for in- 
stance with observations presented in Fig. 3(a) in 
Ref.— This example indicates that in the case of 
a film thickness, which is not an integer multiple 
of Xe, one may observe a spatial variation of the 
number of lamellae in a block copolymer film. 



Figure 21: The spatial coexistence of structures 
with four and five lamellae parallel to the bound- 
aries at a distance Ly. = 4.4373Ae is shown in the 
case of selective boundary conditions i/Zq = V^l, ~ 
1. Parameters £ = g = I and a = 0.015. 

Also the free energy of parallel and perpendic- 
ularly oriented lamellae can be equal for certain 
boundary conditions and film thicknesses, as dis- 
cussed in Sec. IH In Fig. [22l we show an example 
for parameters, where parallel and perpendicularly 
oriented lamellae have the same free energy. Ac- 
cording to the interface between both orientations 
the free energy of the structure in Fig. |22]is slightly 
higher than that of the pure parallel or perpendic- 
ular orientation. Nevertheless, as the interface be- 
tween coexisting lamellae orientations in Fig.l22]is 
not moving the coexisting pattern is long lasting. 




Figure 22: The spatial coexistence of structures 
with lamellae parallel and perpendicular to the 
boundaries at a distance Ly = 4.3854Ae is shown 
in the case of selective boundary conditions i/Zq = 
= 0.4. Parameters £ = g = 1 and a = 0.015. 
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5.3 Dynamics of orientational order- 
ing 

In diblock copolymers confined between bound- 
aries the rotational symmetry is broken. The struc- 
ture factor 5(k, as well as the characteristic 
length scales and ly, as introduced above, do 
not provide a sufficient quantitative characteriza- 
tion of the orientational order of lamellae. Besides 
the so-called Euler characteristics or a complex 
demodulation method— the lamellar morphology 
may be described in the framework of a network 
analysis,-^2jM as applied in this section. 

A basic element of the following analysis is 
image processing and the open source library 
OpenCV— is used for the detection of the inter- 
faces between the A- and 5-rich regions in the 2- 
dimensional binarized images of the field y) . 
The curves along these interfaces are approxi- 
mated by polygonal chains resulting into a set of 
segments of lengths /,■ with the corresponding seg- 
ment orientation angle 0, relative to the x-axis. 
These data allow the calculation of an average in- 
terface segment length 

(^.) = ^I//, (63) 

i=l 

over which neighboring lamellae are parallel to 
each other, as well as the calculation of the aver- 
age orientation of segments and the number of seg- 
ments. These criteria offer an improved distinction 
between patterns of different morphology. 

An order parameter of the segment distribution, 
similar to the order parameter in nematic liquid 
crystals,— is an appropriate quantity for a char- 
acterization of the lamellar patterns during mi- 
crophase separation. Since the orientation angles 
Oi and 0i + TZ are equivalent, the order parameter 
is given by a symmetric second rank and traceless 
tensor 

f) = ( 2xjc Qxy \ 

^ _ Lti/,cos(20,) 

Qxy = z;^— • (64) 



With the scalar order parameter S and the averaged 
orientation angle Q with respect to the x-axis, 

S=^^^^,y. = ^arccos(e,,/5) . (65) 

With the unit vector (the director) n = {nx.riy) = 
(cos 0, sin 0) the tensor order parameter can also 
be written in the following form: 

Qij = S{2nin j - 5ij) , (66) 

where S has for perfectly ordered segments the 
value 5=1 and for an isotropic orientational dis- 
tribution of the segments one has 5 = 0. 

For the same parameters as used in Fig. [Island 
Fig. \T9\ now the quantities 5, and {Q are cal- 
culated as a function of time and the results are 
shown in Fig. [231 The snapshots in Fig. [T8] in- 
dicate a significantly higher orientational order of 
the lamellae in the case of selective boundary con- 
ditions compared to neutral ones. This difference 
in the orientational order can now be quantified by 
comparing S{t) for the two boundary conditions, 
as can be seen in Fig.[23ta). 

The temporal evolution of the average boundary 
segment length (Zy) after a deep quench is shown in 
Fig.[23tc) and the mean orientation of segments 
in Fig. f23\ h) for the three different boundary con- 
ditions: symmetric selective, neutral and mixed. 

As can be already seen in Fig. [TSl the average 
segment length of straight lamellae without de- 
fects takes in the thin film geometry its smallest 
value in the case of neutral boundary conditions. 
Simultaneously, one observes for neutral bound- 
ary conditions the smallest values of the scalar 
order parameter S{t) on the time scale shown in 
Fig. f23\ a.) as well the strongest fluctuations of 
0{t). Since the action length of the boundaries in 
the case of neutral boundary conditions is small 
and the removal of defects is a slow process, (/^) 
increases only slowly as function of time {{Is) 
Ly in the long time limit). 
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Figure 23: In part (a) the scalar order parame- 
ter S, in (b) the lamellae orientation angle 6 [cf. 
Eq. (|65l) 1. and in (c) the averaged segment length 
(Is) [cf. Eq. (|63]) 1 is shown in the strong segre- 
gation regime as a function of time for the same 
parameters as in Fig. [TSland Fig. [T9]for the bound- 
ary conditions symmetric selective (V^ = V^l,. = 
1), symmetric neutral (V^o = Vl,. = 0) and mixed 
(V/o = 0, v/L, = 1). 

In the case of selective boundaries the action 
length of the substrates is larger and therefore ori- 





Figure 24: Snapshots of \j/{x,y,t) at t = 100 (a), 
t = 200 (b) and t = 300 (c) are shown for sym- 
metric selective boundary conditions and the same 
parameters as in Fig. \TS\ 

ented lamellae are formed much earlier. Conse- 
quently, one observes higher values of (Is) and 
of the order parameter S much earlier. The fact, 
that S{t) has still not reached the value 5 = 1 at 
about t — 10^ in Fig.[23la) is related to the few de- 
fects left, as can be seen for instance in Fig. [TST b). 
As the regular structure is represented by lamel- 
lae parallel to the substrates, (h) — )■ Lx for long 
time dynamics. The reduction of (Is) for selec- 
tive boundaries at about t = 200 in Fig. f23tc) is 
related to intermediate structures that occur during 
coarsening as indicated by the transition from the 
pattem in Fig. l24T a) to the pattern in Fig. |24tb). 

The rather early achieved orientational order for 
selective boundaries is also indicated by the be- 
havior of 0{t), which approaches zero quite early 
in Fig. [23lfb). The large fluctuations of 6{t) in 
Fig. \23\ h) for the neutral boundary conditions re- 
flect the coarsening and the related removal of de- 
fects on the route to a higher orientational order. In 
this case the boundary segments are preferentially 
oriented perpendicular to the substrates resulting 
into the orientation angle ^ 7r/2 for long-time 
evolution. In case of mixed boundaries, one ob- 
tains a mixing of both trends with respect to the 
orientation of lamellae. The selective (upper) sur- 
face triggers lamellae oriented parallel to the sub- 
strate whereas the neutral (lower) surface initiates 
lamellae oriented perpendicular to the substrate 
with less defects as for the two neutral boundaries 
[see Fig. [T8lb)1. and accordingly the results for S 
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and (/v), represented by dotted lines in Fig. U3\ a) 
and Fig. [23lc). lie between the two symmetric 
cases. 

For comparison we show the temporal evolution 
of S, 6 and (Is) in Fig. [25] in the weak segrega- 
tion regime with £ = 0.37. Both, the behavior of 
S{t) and of (Is) confirm, that in the weak segrega- 
tion regime an orientational order is reached on a 
smaller time scale as in the case of the strong seg- 
regation regime. 

The composition waves with almost equilibrium 
wavelength result into {Is) ~ Lx for the selective 
and mixed boundary conditions. The selective 
boundaries provide the fastest formation of regu- 
lar parallel lamellae and therefore the fastest sat- 
uration of (Is). In case of neutral boundaries {Is) 
is increased faster in time compared to the deep 
quench [Fig.[23tc)] with the tendency (Is) — )■ Ly. 

In case of a not too deep quench the scalar or- 
der parameter S shown in Fig. [25l a) grows much 
faster in time compared to the deep quench. The 
boundary segments are again preferentially ori- 
ented perpendicular to the boundaries for the neu- 
tral boundary conditions and parallel to the bound- 
aries for the selective and mixed boundary condi- 
tions [Fig.|25tb)]. 

Thus even this quite simple analysis of lamel- 
lar patterns provides quantitative characteristics of 
the dynamics and influence of the boundary con- 
ditions on the resulting patterns, that are compli- 
mentary to the standard analysis of the structure 
factor. 
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Figure 25: In part (a) the scalar order parameter 
S and in (b) the lamellae orientation angle Q [cf. 
Eq. (l65l)1. and in (c) the averaged segment length 
(/j) [cf. Eq. (I63l) 1 is shown in the weak segregation 
regime as a function of time for the same parame- 
ters as in Fig.|20]for the boundary conditions, sym- 
metric selective (v/q = V^l,. = 1)^ symmetric neutral 
(V^ = V^L, = 0) and mixed (v/q = 0, \\fL, = 1). 
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6 Summary and conclusions 

The formation and stability of lamellae in block 
copolymers has been investigated in terms of a 
mean-field model. A method for the determina- 
tion of the stable wave number band has been in- 
troduced in Sec. |3] and the shape of this band pro- 
vides the basis for a deeper understanding of stable 
lamellae conformations which are locked in ex- 
periments on BCP films to different wavelengths 
when using spatially periodic chemical nano pat- 
terning of substrates.— 

We have found, similar as in previous calcula- 
tions in terms of self-consistent mean field theo- 
ries or phenomenological free energy models,-^^*^ 
that selective boundaries induce lamellae orienta- 
tions parallel to the substrates and in the case of 
neutral boundaries lamellae orient perpendicular 
to surfaces. We present also estimates in terms of 
the different length scales parallel and perpendicu- 
lar to the lamellae, whether lamellae orient parallel 
or perpendicular to confining boundaries. Some of 
the lamellae conformations calculated within this 
work resemble very much the lamellae orienta- 
tions observed experimentally in thin BCP films at 
neutral substrates in Refs.— We derive also ana- 
lytical expressions in the case of lamellae parallel 
to substrates for the concentration modulation per- 
pendicular to the boundaries, which can be useful 
for qualitative considerations in further works. 

In the case of mixed boundary conditions, i.e. 
selective at one boundary and neutral at the oppo- 
site boundary, we find a critical value \ifs{crit) of 
the selectivity below which the energetically pre- 
ferred, homogeneous lamellae orientation changes 
from parallel to perpendicular with respect to the 
confining boundaries for any film thickness. 

The results obtained are interesting also with re- 
gard to a recently used strategy to control the long 
range lamellae order in diblock copolymer films, 
where a thickness-dependent orientation of lamel- 
lae has been found.— While the lamellae oriented 
parallel to the substrate in the ranges of film thick- 
nesses, J < 19 and J > 23, a perpendicular orienta- 
tion in the range 19 < J < 22 was observed.— At a 
first sight, this experimental observation seems to 
be in contradiction to our results. However, these 
experiments were done in the presence of a sol- 
vent that changed the degree of swelling of the 



BCP film. Although our mean-field model does 
not contain explicitly the effects of a solvent on 
the lamellae formation and its interaction with sur- 
faces, we suggest to use our results for an interpre- 
tation of the mentioned experiment. It has been 
found that in the presence of a solvent the degree 
of swelling ^ = Xe/K (where Xs stands for the 
lamellae period in the swollen state) depends on 
the film thickness. For very thin films {d ~ 3) the 
degree of swelling is around ^ 0.68 whereas for 
thicker films {d ~ 30) it is about ^ 0.715, which 
means, thicker films swell about 4% less than thin- 
ner films and the concentration of the solvent 
is decreasing with increasing the film thickness. 
Accordingly, for thinner films in the swollen state 
we calculate an "effective value" of the model pa- 
rameter a ~ A^'^ : as = (0.68)^ao (cto indicates 
the interaction parameter in the absence of a sol- 
vent). It is about 18% smaller than for thick films 
where Us = (0.71 5)"^ ao- Thus to model the influ- 
ence of the solvent we can assume that the param- 
eter a is effectively increased by increasing the 
film thickness. In addition our simulations reveal 
that the critical absorption at the surface, \ifs{crit), 
increases with a. Assuming a linear behavior of 
\lfs{crit) as a function of a 

\\fs{crit) = a-a + b, (67) 

in a small range around as{d = 30) /as{d = 3) ^ 
1.2, we calculate the slope a ^ 14.8. Therefore, in 
the range J = 3 • ■ ■ 30, 

y/s{crit) [d = 30] - Ysicrit) [d = 3]^ 0.7ao , 

(68) 

and \ifs{crit) increases slightly with the film thick- 
ness. 

On the other hand the solvent may reduce the 
selectivity of the confining surface (see, e.g.,—). 
Thus, is decreasing with an increasing solvent 
concentration. According to the swelling behav- 
ior it means that i//^ is increasing with increasing 
the film thickness d as in thicker films the solvent 
concentration is lower. The exact form of the curve 
^fs{d) can be determined experimentally by mea- 
suring the BCP-substrate interfacial tension in the 
swollen state for various film thicknesses. 

Combining now the two effects, the addition of 
a solvent has, we end up with two curves that, de- 
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pending on their relative position, offer the possi- 
bility for a reorientation effect of the lamellae as a 
function of the film thickness (see Fig.[26l). The re- 
gion where x^^sid) < Ws{crit) indicates a possible 
transition region. 
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Figure 26: Critical value y/sicrit) of the surface 
selectivity and i//^ rescaled by CUq as a function of 
d. The additional axis at the bottom indicates the 
solvent concentration c.,. 

Of course the curve for \]/s{d) in Fig. |26]is to 
some extend hypothetical and the shape has to be 
determined experimentally via measuring the in- 
terfacial tension. Nevertheless, our results indicate 
a route for further experiments on the thickness- 
dependent lamellae reorientation. Especially for 
diblock copolymers with a pronounced change of 
the degree of swelling as a function of the film 
thickness (like, e.g., in,— where the increase of 
the solvent uptake with decreasing film thickness 
is more than 10%) a considerably wide reorienta- 
tion range may be realized by a suitable tuning of 
the wetting properties of the confining surfaces. 

In addition to the energetic considerations of 
the influence of boundaries on the homogeneous 
lamellae orientations, we also investigated the dy- 
namical evolution of lamellae structures between 
boundaries. In the case of mixed boundaries, one 
also finds complex lamellae conformations, even 
if they have a higher free energy than a homoge- 
nous lamellar order either parallel or perpendic- 
ular to the confining parallel boundaries. Simu- 
lations of the time-dependent mean field model 
show, that the type of boundary condition deter- 
mines strongly the evolution of the orientational 



order as well as the number of defects in BCP 
films parallel and perpendicular to the boundaries, 
which has been quantified by using an order pa- 
rameter for the characterization of the lamellae ori- 
entation. The consideration of different quench 
depths in combination with various boundary con- 
ditions provides a strategy for the experimental 
preparation of defect-free oriented lamellae. 
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A Numerics of nonlinear solu- 
tions 

Here we describe the numerical determination of 
stationary periodic solutions of Eq. (fT2] | for /3 = 
as well as their linear stability. Since Eq. (fT2] ) is 
isotropic one can choose the x direction parallel to 
the wave vector of the periodic solution. A Fourier 
expansion of the periodic solution y^ki^), as given 
by Eq. (|39l ). leads together with Eq. (fT2] ) after pro- 
jection onto the j-th Fourier mode e'*^-' to a set of 
nonlinear equations for the coefficients Af. 

[£{jkf - {jkf - a\Aj - {jkfY.^iA,„Aj-i-'n = , 

l,m 

i = -M...M. (69) 

For M > 1 this system of nonlinear equations is 
solved numerically by Newton's iteration method 
and M is adjusted to keep the relative error smaller 
than 10^. For £ = 1 this accuracy can be main- 
tained in the case of a = 0.001 with a very steep 
density profile by choosing M = 120 modes and in 
the case of smoother density variations by M = 15 
modes. For larger values of a a smaller number 
of modes is required as the solution becomes more 
harmonic. 

The linear stability of periodic solutions \ifk{x) 
of Eq. (fT2l) with respect to small perturbations 
i//i(jc,j, f) is investigated as follows. One starts 
with the ansatz i//(a:) = Wk{^) + Wii^iJ-:^) ^i^d a 
linearization of the basic equation (fT2] ) with re- 
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spect to v/i t) gives 

d,\i/i {x,y,t) = V2 {-e + Sxj/^ - V-) xj/i - axj/i ,(70) 

wherein the spatially periodic function y/k enters 
parametrically. For a solution of this linear equa- 
tion (TTOl) with periodic coefficients one uses a Flo- 
quet ansatz, 

xi/i{x,y,t) = e^^e"(^'^°^^+>'^i"^)<^^(x) +c.c. , (71) 

with the Floquet parameter s, a 2;r//c-periodic 
function (^f (x) and the angle enclosed between 
the wave vector of the perturbation and the 
wave vector of the basic periodic solution. This 
periodic function ^f{x) can be represented by a 
Fourier expansion 



M 



ikxn 



(72) 



n=—M 



Taking into account Eq. (|39| ) for x^fk the linear par- 
tial differential equation Eq. (iTOl) is transformed 
after projection into an eigenvalue problem 



oDn = {eCn- Cl -a]Dn- 3C„ Y^MA^D^ 



Cn = {kn + s COS ey + s^sm^e 
n = -M .. .M , 



(73) 



where the coefficients A/ are determined by 
Eq. ( |69| ). We are interested in the growth rate 
o{e,k,s,Q), i.e., in the eigenvalue o with the 
largest real part. The condition Re[o{e,k,s, 9)] = 
yields the stability boundaries £ = £o(fc, 5,0). 
For Re{o) = and 0=0 one finds £ = £E{k) 
that determines the Eckhaus boundary. In the case 
of Re{o) = and = k/2 the corresponding 
£ = £zz{k) gives the zig-zag line. 



Stability of weakly nonlinear 
solutions 



B 



As described in Sec. 13.41 a periodic solution in a 
two-dimensional isotropic system may be destabi- 
lized by modulations along the wave vector (Eck- 
haus instability), undulations perpendicular to it 
(zig-zag instability), and a combination of both 



types (skewed varicose).— Two of the instability 
branches are given in Fig. [Hand they may be de- 
termined analytically near threshold by analyzing 
the stability of the solution given by Eq. (l35l) for 
P = with respect to small perturbations 5 A: 



(74) 



The analytical form of the perturbation is as fol- 
lows: 

dA = e^'e'^^'iaiexplisi^xcosG +ysmG)] 

-|-a2exp[— Z5(xcos0 +ysin0)]} . (75) 



A linearization of Eq. (1331 ) with respect to small 
perturbations 5A leads for the growth rate a to 
a 2 X 2 eigenvalue problem (see e.g.^^'^^). In the 
special case = and under the condition of neu- 
tral growth a = the following expression for 
the control parameter at the Eckhaus stability 
boundary (0=0, longitudinal instability) follows: 



6^ 

kl 



6{k-r 



(76) 



A comparison with the neutral curve in Eq. (I37l 
shows that the width of curve rE{Q) is narrower 
than rt^{Q) by the famous factor 



1 



rE{Q) 
rN{Q) V3 



(77) 



for Eckhaus stability boundary .-^Sj^^d^iH 

The zig-zag stability boundary (transversal in- 
stability) results for the case = n/2 and a = 0: 



2 = 0, i.e., fez = 1 . 



(78) 



Hence, the stationary weakly nonlinear solution 
given by Eq. (l35l) is linearly stable in the region 
Q>0(k> 1) between the zig-zag line in Eq. (1781) 
and the Eckhaus boundary Qe = a/ rk^ / 6 . 

The results for the stability boundaries found so 
far in this appendix in the framework of the ampli- 
tude equation Eq. (|33l ) are typically valid only in 
the vicinity of the critical point (r > 0, |fc — fe| ^ 
1). An essential improvement of the analytical re- 
sults for the stability diagram can be obtained by 
the Galerkin approach in a one-mode approxima- 
tion, as described in the following. 
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Inserting the ansatz given by Eq. (|32l ) into 
Eq. (|3T1) and projecting on the critical mode 
one gets in the leading order 

d,A = Vl\-2kl{r+ i)A + 3\A\^A-VlA]-klA , 



with 



k. 



(79) 



Above threshold one may derive for the amplitude 
Aq with the ansatz Eq. (I35l) via Eq. ( 1791) the follow- 
ing expression: 



2l2- 



[{Q + k,)^- k^] 
2k^.iQ + kc 

(P-2)2^ 
2P 



(80) 



One can easily see that this amplitude Aq van- 
ishes at the neutral curve Eq. (|28T) for arbitrary val- 
ues of r. Inserting perturbation Eq. (1741 ) with the 
amplitude given by (|80l ) into Eq. (1791 ) the growth 
rate o is again calculated from a 2 x 2 eigen- 
value problem. The various stability boundaries 
are determined via the neutral stability condition 
a(r, 2,0) = in terms of the control parameter 
r(2, 0) by keeping simultaneously only the lead- 
ing terms in s. Minimization of r{Q, 6) with re- 
spect to 6 gives in the range Q>0 the angle 0=0 
and therefore an Eckhaus stability boundary 



rE 



3P - 2k^ -6P + 5 
2^(^4 + 3) 



(81) 



that coincides with the result obtained in the 
framework of the free energy considerations [see 
Eq. (l49l)1. For k between the neutral curve given by 
Eq. (l28l ) and the Eckhaus boundary in Eq. (ISTT ) the 
periodic solution in Eq. (l35T l with Aq from Eq. (l80l l 
is unstable with respect to long-wavelength per- 
turbations along the wave vector, i.e. s ^ and 
0=0. 

In the vicinity of the band center k = 1 one has in 

the leading order rE = 6(k—l)^-\ in agreement 

with the result derived via the standard amplitude 
equation [see Eq. (1761)1. 

A minimization of r(2, 0) in the range Q < 
gives the angle = ;r/2 for the zig-zag instability 
line 



in agreement with Eq. (1501) . For the perturbations 
Eq. dill) with s ^ 0, 6 = 7t/2the growth rate o is 
negative for k on the right hand side of the zig-zag 
line Eq. (l82l) up to the Eckhaus boundary Eq. (ISTl) . 

A skew varicose instability with < < ;r/2 
does not occur. 



C Weakly nonlinear solution 
under confinement 

Close to the onset of microphase separation (r > 
0) the dynamics of the amplitude of the periodic 
order parameter field \j/{r) is govemed by the 
Newell-Whitehead-Segel amplitude equation (l33l) . 
as described in Sec. 13.21 This equation has spa- 
tially periodic solutions in extended systems as de- 
scribed in Sec. 13.21 but it may also be solved in the 
presence of boundaries. 

Here we take into account boundary conditions 
for the case of lamellae oriented parallel to the sub- 
strates. With the ansatz 



^5(3;)e*^(>--^'')+c.c.. 



(83) 



one gets, starting from Eq. (1311 . the following 
equation of the envelope B{y): 

^d^B+{l-\B\^)B = 0. (84) 

This equation has constant solutions of the form 

5o(y)=5o = ±l, (85) 

corresponding to a spatially periodic field '^f{y) of 
constant amplitude. 
Eq. ( l84l ) also has the solution 



5(>;) =5otanh 



iy-y) 



(86) 



and this may be used to construct approximate so- 
lutions for lamellae parallel to the two boundaries. 
The decomposition (l83l) is based on the assump- 
tion, that B{y) varies slowly on the scale 271 /kc and 
close to threshold one may simplify the boundary 



Qzz = , i.e., kzz = 1 



(82) 
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conditions ( |59l ) to the following conditions, 

dyW\y=o,Ly = 0, (87a) 

W\y=o = ¥o = c-V/b, (87b) 
V/|v=L,, = V^L, = c- V/fc, (87c) 

by neglecting higher order derivatives of i//. The 
constant c is used to switch between different types 
of boundary conditions: c = corresponds to neu- 
tral boundaries, and c ^ to selective (symmetric) 
boundaries. 

In order to fulfill the boundary conditions (ISTT l at 
y = and y = Ly, we use a linear superposition of 
the solution (l86l) of the form 



(88) 



B{y)=vi+V2Yi{y)+v^Y2{y). 



where the constants vi,V2,V3 indicate the possible 
signs ± 1 and with 



Yi {y) = tanh 
F2(y) =tanh 



iy-yi] 



kc\/r 



(y-yi) 



(89a) 
(89b) 



The constants yi and y2 are determined by the 
boundary conditions. 
With the bulk value 



Skh 



(90) 



the order parameter \i^{y) takes the form 



¥{y) = Vb[vi+V2Yi{y) +V3Y2{y)]cos[kc{y -yo)] . 

(91) 

Inserting Eq. (|9TI) into the boundary conditions 
(f87l) one has the following equations 



V2(l-r2(0))^ 



cos{kcyo)+ 
vi + V2Y1 (0) - V3] sm{kcyo) = , 

cos(fccjo) + 



V3(l-y2n^y))^ 



[vi +V2 + V3y2(^)]sin(/:cyo) =0, 
[vi +V2Yi{0) -V3]cos{kcyo) =c, 

[Vi -|-V2+V3y2(^)]cos(^c3'o) =C 



(92a) 



(92b) 
(92c) 
(92d) 



Here we assumed a system size Ly = n ■ {2n)/kc 
and exploited the approximation that the bound- 
aries do not influence each other so that Yi (Ly) — 1 
and ^2(0) = —1. From this one may directly de- 
duce Vi = V3 = — V2. 

In the following some special cases are con- 
sidered to illustrate the high quality of this ap- 
proximation. The cases of selective (c = 1) or 
neutral (c = 0) boundaries may be solved explic- 
itly. For the case of neutral boundaries we obtain 
kc-yo = ±7t/2, 3^1 = and y2 = Ly as a solution [see 
Fig. |9](b)]. vi remains arbitrary in this case. If 
c = 1 the value at the boundary corresponds to the 
bulk value and this case is conform with the pe- 
riodic one where B{y) = ±1 = const, and yo — 
[see Fig. |27](a)]. For c < 1 the selectivity is re- 
duced [see Fig. [27] (b)]. In this case an explicit so- 
lution is not possible but one may use the approx- 
imation s'm{kcyo) ~ if c is not too small what 
corresponds to yo ~ 0. With this approximation 
we find ^1(0) ^ c and Y2{Ly) ^ — c and vi = — 1 
as a solution. The example shown in Fig. [27] (b) 
is for c = 0.5 what results in yi f» —10.83 and 
y2 ~ Ly + 10.83. Although this is only an approx- 
imation it still fits very well the full numerical so- 
lution. Of course this case may be studied more 
accurately by taking into account the wave num- 
ber change. In this case one has more free param- 
eters to adjust reasonably. Furthermore in the case 
of mixed boundaries this approximation becomes 
more complex and higher order derivatives have to 
be taken into account. The solution of the result- 
ing equations becomes involved and the advantage 
of the approximations gets lost in comparison with 
the full numerical solution of the problem. 



for the determination of the constants of the ansatz. 
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Figure 27: The numerical solution of 
Eq. (fT2l) is shown for various boundary conditions 
for lamellae parallel to the boundary: In part (a) 
for selective boundaries c = 1 and in (b) for selec- 
tive boundaries with a reduced selectivity c = 0.5. 
The solid lines mark the numerical solution and 
the analytical approximation given by Eq. (1911 is 
displayed by the dots. The parameters are Ly = 
lOAg, g= l,r — 0.021, kc = 0.7 (corresponding to 
£ = 1 and a = 0.24). 
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